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Abstract. We present differential radial velocity mea- 
surements of Proxima Centauri collected over 4 years with 
the ESO CES with a mean precision of 54 ms" 1 . We find 
no evidence of a periodic signal that could corroborate the 
existence of a sub-stellar companion. We put upper limits 
(97% confidence) to the companion mass ranging from 1.1 
to 22 Mj up at orbital periods of 0.75 to 3000 d, i.e. sep- 
arations 0.008 — 2 AU from Prox Cen. Our mass limits 
concur with limits found by precise astrometry (Benedict 
et al. 1998a and priv. comm.) which strongly constrain the 
period range 50 — 600 d to 1.1 — 0.22 Mj up . Combining 
both results we exclude a brown dwarf or supermassive 
planet at separations 0.008 — 0.69 AU from Prox Cen. We 
also find that, at the level of our precision, the RV data 
are not affected by stellar activity. 
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1. Introduction 

Searches for companion objects to Proxima Centauri 
(Prox Cen, GJ551, V645 Cen, HIP70890), the nearest star 
(d = 1.2948 ± 0.041 pc; M5Ve), date back 20 years to the 
astrometric work by Kamper & Wesselink (1978) and the 
infrared photometric scanning by Jameson et al. (1983). 
The most extensive search so far consists of astrometric 
monitoring with the HST FGS #3 (Benedict et al. 1998a) 
where an astrometric precision of 0.002" per axis and 
a detection limit for an astrometric variation of 0.001" 
was achieved, thereby strongly constraining the mass of 
a possible companion. Limits to the iVband magnitude 
of objects within projected separations of 1 — 10 AU 
from Prox Cen were found by Leinert et al. (1997) to be 
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K = 12.9 — 15.1 mag, i.e. AK = 3 — 5 mag below the 
empirical end of the main-sequence. 

Recently the efforts to detect substellar objects near 
Prox Cen culminated in the announcement of a possible 
companion by Schultz et al. (1998) who used the HST FOS 
as a coronographic camera. These authors reported excess 
light near Prox Cen seen in two images separated by 103 d. 
Within this time span, the suspected object appeared to 
have moved in separation from 0.23" to 0.34" (indicating 
a separation near 0.5 AU) and in P. A. from 45° to 100°. If 
interpreted as a companion the object would be ~ 7 mag 
fainter in the FOS red detector. 

However, a subsequent observation with the HST 
WFPC2 at two epochs (separated by 21 d) by Golimowski 
& Schroeder (1998) could not verify the existence of any 
companion object to Prox Cen within a separation from 
0.09" to 0.85" (0.11 - 1.1 AU). The authors concluded 
that they should have seen the object at only sa 3.7 mag 
fainter than Prox Cen in their images taken at l^tm. Con- 
sequently, they suggested that the excess light seen by 
Schultz et al. (1998) was an instrumental effect. 

In this paper we report on 4 years of precise radial 
velocity (RV) monitoring of Prox Cen and contribute new 
evidence to the debate on a substellar companion. 

2. The planet search program at the ESO CES 

Our planet search program of 39 late-type stars us- 
ing high-precision RVs was begun at ESO La Silla in 
Nov. 1992. Prox Cen was first observed in July 1993. We 
used the ESO 1.4m CAT telescope and CES spectrograph 
equipped with the f/4.7 Long Camera and ESO CCDs #30 
or #34. The obtained resolving power, central wavelength 
and spectrum length were 100, 000, 5389 A, and 48 A. 
For high measurement precision for differential RVs we 
self-calibrated the spectrograph using an iodine (I2) gas 
absorption cell temperature controlled at 50° C (Kurster 
et al. 1994; Hatzes et al. 1996; Hatzes & Kurster 1994). To 
obtain RV measurements we model the stellar spectra as 
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observed through the iodine cell using a 'pure' stellar spec- 
trum (recorded without the iodine cell in the light path) 
and a 'pure' iodine (7 2 ) spectrum from dome flat measure- 
ments. The resulting RV data are then corrected to the 
solar system barycenter via the JPL cphcmeris DE200. 

For stars brighter than 5.5 mag our short-term (i.e. sin- 
gle night), best case long-term, and working long-term 
precisions (i.e. obtained under all observing conditions) 
are 4—7, 11, and 20 — 25 ms^ 1 , respectively (Kiirster 
et al. 1994, 1998, 1999). Prox Cen is unique in our sam- 
ple in that it is by far the faintest star we have observed 
(V = 11.01 mag). 

3. The RV measurements for Proxima Cen 

Tab. 1 shows the journal of observations of our differential 
RV measurements. A total of 58 spectra from 29 nights 
were available. Before further analysis the spectra were 
combined into nightly bins (col. 1) as outlined below, with 
the bins containing between 1 and 5 spectra (col. 2). 

While self-calibration with an iodine cell is an excel- 
lent method to overcome instrumental instabilities such 
as instrumental drifts other instabilities such as focus and 
alignment changes or instrument vibrations require (in 
principle) additional modelling of the instrumental pro- 
file (IP; Butler et al. 1996; Valenti et al. 1995). At low 
signal-to-noise (S /N) ratios such as obtained for Prox Cen 
IP reconstruction becomes impossible. However, to some 
extent one can overcome the remaining measurement un- 
certainty by modelling the observed spectra (star+iodinc) 
with various combinations of pure star and iodine spectra. 

Ten different pure star and pure iodine pairs (of the 
same night) were available to build models for each of the 
58 star+iodine spectra. Based on goodness-of-fit some of 
the most inadequate models were rejected. Col. 3 gives 
the total number of the accepted star+iodine models for 
all the spectra in each nightly bin (sum over all spec- 
tra and models). For an intercomparison of the different 
models their RV zero points have to be matched. To do 
this we subtracted for each model the mean RV for the 
whole times series. Subsequently, we averaged for each 
star+iodine spectrum the RV data from the individual 
models. At last, nightly (bin) averages for the Julian day 
(col. 4) and RV data corrected to zero mean (col. 5) were 
calculated. 

A first estimate of the RV error (col. 6) was then based 
on the rms scatter in each data bin together with a prop- 
agation of the error introduced by the process of match- 
ing the zero points for the different models. However, we 
found a positive correlation with a correlation coefficient 
of +0.514 between the total number of models in a bin 
(col. 3) and the error estimate (col. 6) meaning that these 
errors tend to be smaller when estimated from smaller 
numbers of models. This indicates that these error esti- 
mates are not representative of the true errors, in partic- 
ular those from small numbers of models. 



Table 1. RV data for Prox Cen. Col. 1: bin number; 
col. 2: number of spectra in the bin; col. 3: total number 
of models for all spectra in the bin (i.e. the sum over all 
spectra and models); col. 4: mean Julian day of the bin; 
col. 5: mean differential RV of the bin; col. 6: rms scatter 
within the bin; col. 7: combined S/N ratio of the bin; col. 8: 



equivalent RV error as estimated from the S/N. 


# # E 


JD - 


ARV 


rms 


comb 


err oq 


Sp 


2, 400, 000 


[m/s] 


[m/s] 


S/N 


[m/s] 


1 1 7 


49179.6643 


106.35 


44.85 


21.46 


71.78 


2 2 14 


49246.5034 


- 18.00 


78.30 


17.88 


86.16 


3 2 18 


49412.8484 


3.48 


46.24 


23.87 


64.54 


4 4 34 


49413.7676 


- 78.83 


68.73 


34.10 


45.18 


5 5 40 


49492.7022 


52.55 


85.96 


39.78 


38.73 


6 4 30 


49548.5800 


56.66 


65.18 


35.12 


43.86 


7 2 15 


49602.5149 


- 2.23 


60.60 


28.44 


54.17 


8 2 17 


49794.7186 


107.13 


103.27 


26.28 


58.62 


9 2 14 


49795.6955 


28.33 


44.09 


28.10 


54.82 


10 2 16 


49906.5524 


0.64 


44.11 


27.91 


55.20 


11 1 9 


49907.5940 


60.47 


38.37 


21.47 


71.75 


12 1 8 


50146.7675 


- 70.90 


42.87 


25.89 


59.50 


13 2 17 


50229.6486 


- 40.73 


56.75 


36.54 


42.16 


14 1 8 


50230.6858 


55.73 


49.09 


26.82 


57.44 


15 2 14 


50304.4963 


2.19 


52.60 


30.92 


49.82 


16 2 15 


50313.5397 


- 35.71 


43.80 


40.39 


38.14 


17 2 14 


50314.5965 


- 37.08 


49.38 


37.62 


40.95 


18 2 13 


50319.5519 


- 21.17 


45.65 


31.42 


49.03 


19 1 6 


50323.5606 


- 88.32 


30.21 


25.87 


59.55 


20 1 7 


50325.5852 


- 47.33 


36.42 


27.10 


56.85 


21 2 17 


50348.5179 


6.32 


45.15 


38.82 


39.68 


22 2 18 


50349.5142 


- 18.19 


49.07 


36.90 


41.75 


23 1 7 


50359.4920 


- 36.71 


42.60 


24.56 


62.72 


24 2 17 


50477.8459 


45.31 


73.70 


25.24 


61.03 


25 2 14 


50478.8507 


-104.59 


93.69 


21.67 


71.09 


26 2 18 


50524.7935 


43.71 


48.86 


28.80 


53.49 


27 2 17 


50551.6664 


44.23 


48.27 


28.03 


54.96 


28 2 12 


50570.7647 


- 3.27 


61.04 


38.30 


40.22 


29 2 15 


50648.5032 


- 10.02 


42.17 


38.09 


40.44 



As an independent estimate of the weight of each data 
bin, col. 7 shows the combined S /N ratio per spectral pixel 
(mean values) of each data bin. We know from simulations 
(Hatzes & Cochran 1992) that the RV measurement error 
ARV oc (S/N)^ 1 which serves to estimate the RV error. 
Choosing the constant of proportionality such that the 
mean of the resulting errors is equal to the total scatter in 
the RV measurements (53.9 ms" 1 ) we obtain the equiv- 
alent RV errors listed in col. 8 of Tab. 1. Fig. 1 displays 
a time series of our RV data (29 bins) for Prox Cen with 
error bars corresponding to these equivalent errors. 

4. Period search 

To look for a periodic signal that could manifest the pres- 
ence of a companion, we searched the frequency range 
/ min = 1/T ... / max = l/(2At), or 0.0007 ... 0.5494 d"\ 
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Fig. 1. Time series of the differential RV data of Prox Cen. 

where T is the total time baseline and At is the minimum 
separation between data points. Thus a period range of 
1.8202 ... 1428.6 d was searched. The choice of the maxi- 
mum frequency was made in analogy to the Nyquist cri- 
terion which, however, is well defined only for equidistant 
data sampling. Since most of our data sampling is much 
cruder than our minimum sampling any signals at periods 
shorter than about 5 d should be treated with care. 

Two different types of periodogram were used: a) the 
Scargle periodogram (Scargle 1982) which is equivalent to 
least squares sine fitting with equal weight for all data 
points; the power is proportional to the square of the am- 
plitude of the corresponding sine wave; b) a sine-fitting 
routine that takes into account data errors by minimizing 
X 2 ; we used the equivalent errors (col. 8 of Tab. 1). 

Fig. 2 shows only the Scargle periodogram, since the \ 2 
minimization approach yielded a very similar result that 
does not change the interpretation. ^From 10,000 runs of a 
bootstrap randomization scheme (see Kiirster et al. 1996; 
Murdoch et al. 1993) we determined the levels of the false 
alarm probability <f> (FAP) corresponding to various power 
levels. As shown in Fig. 2 all periodogram peaks are in- 
significant having > 90%. 

In contrast to searching a period range signals at a 
priori known periods can be significant at smaller power 
levels z. For the Scargle periodogram the FAP is given by 
<p = 1 — (1 — e~ z )™, where n is the number of independent 
frequencies in the search interval (Scargle 1982). Hence for 
a signal at a single a priori known period <ft — e ~ z ■ 

Periods that may be present are the period of the stel- 
lar rotation and that of the activity cycle. RV searches 
strongly benefit from ancilliary information such as the 
knowledge of these periods that aids the interpretation 
of RV data. Star spots and inhomogeneous granulation 
patterns in active stars cause distortions in stellar absorp- 
tion lines; when sufficient resolution and/or signal-to- noise 
is lacking these distortions can be mis-interpreted as RV 
shifts that vary with the rotation period (rotational modu- 
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Fig. 2. Periodogram for the RV data. Horizontal dashed 
lines correspond to various levels of the false alarm prob- 
ability. 

lation) or with the activiy cycle. Estimates for the rotation 
period of Prox Cen were given by Benedict et al. (1998b; 
based on HST FGS photometry) finding P Iot — 83.5 d plus 
variability at the first harmonic, and by (Guinan & Mor- 
gan 1996; monitoring of the Mgll h+k flux with IUE) who 
found 31.5 ± 1.5 d. Benedict et al. (1998b) also estimate 
the length of Prox Cen's activity cycle to be « 1100 d. 

We do not find significant FAPs at any one of these 
periods nor at first harmonics or periods twice as large 
(relevant in case the literature values are first harmon- 
ics themselves). Allowing for some uncertainty in these 
period values we also searched their vicinities finding 4> = 
0.47% (0.16%) for P = 150.4 d and <j> = 1.32% (0.92%) for 
P = 30.4 d, where theoretical values, <j> = e~ z , as well as 
values bootstrapped from 10,000 runs (values in brackets) 
are given. At best, i.e. only if one allows for considerable 
errors in the original period estimates, a marginal detec- 
tion of (1) a period twice as long as the rotation period 
by Benedict et al. (1998b) and (2) the rotation period by 
(Guinan & Morgan 1996) might be indicated. 

Reconfiguration of active regions causes amplitude and 
phase changes in rotationally modulated signals compli- 
cating their detection in data that extend over many ro- 
tation cycles. However, our results concur with the find- 
ings by Saar et al. (1998) that predict an activity-related 
RV scatter of < 10 ms _1 for rotation periods > 16 d. It 
appears that the RV variation seen for Prox Cen is repre- 
sentative of our measurement precision for this faint star, 
and cannot be attributed to instrinsic stellar variability. 

5. Limits to companion parameters 

Lacking a clear RV signal in our Prox Cen data we used 
a Monte Carlo simulation to derive upper limits to the 
mass of still possible companions in the period range 
0.75 — 3000 d. Random data sets with the same tempo- 
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Fig. 3. RV derived upper limits to the mass of companions 
to Prox Cen as a function of orbital period and separation. 
The four parallel curves represent the to sin i and several 
confidence intervals for the mass limit. The single bold 
line is the mass upper limit from HST FGS astrometry 
(Benedict, priv. comm.). 



ral sampling and with the rms scatter of our data were 
created. We added sinusoidal signals with different pe- 
riods, amplitudes, and phases, and evaluated their peri- 
odograms. At each period the amplitude was determined 
for which 99% of the periodograms showed a power cor- 
responding to > 1% FAP. Hence the combined confidence 
for the detection of a sinusoidal signal is 98%. For eccentric 
orbits it may be somewhat lower. 

Fig. 3 shows the to sin i (to the planet mass, i the or- 
bital inclination), that corresponds to RV amplitudes at 
this confidence level, as a function of period and sepa- 
ration. We assumed a stellar mass of 0.12 M Q (for an 
M5Ve star; Kirkpatrick & McCarthy 1994). To account 
for the unknown inclination one can use its probabil- 
ity distribution (for random orientation of the orbits). 
The probability that i exceeds some angle 9 is given by 
p(i > 9) = cos(#). From this one can construct confidence 
intervals for the true companion mass. There is a confi- 
dence of 90% (95%, 99%) that the true mass is no more 
than a factor 2.294 (3.203, 7.088) larger than the to sin i. 
Mass limits for these confidence intervals are also included 
in Fig. 3 with the corresponding labels. Values in brack- 
ets are the combined confidence levels accounting for both 
the confidence of the inclination and the confidence of de- 
tection. We choose the curve with the highest confidence 
(97% combined) as the RV derived upper mass limit. 

We also got upper limits (99% confidence) to the com- 
panion mass from HST FGS astrometry (Benedict, priv. 
comm.; cf. Benedict et al. 1998a) that are included in 
Fig. 3. Being most stringent at longer periods they comple- 
ment our RV derived limits constraining the period range 
50 — 600 d. When combining them with our RV mass lim- 



its we can exclude massive planets around Prox Cen over 
a wide period range as detailed in Sect. 6. 

6. Conclusions 

1. Prox Cen does not have a close (w 0.4 AU) brown 
dwarf companion as suggested by Schultz et al. (1998). 

2. RV derived upper mass limits range from 1.1 to 
3.7 Mj up for periods from a few days to a few weeks. 

3. In the period range 50 — 600 d (separations 0.13 — 
0.69 AU) the RV derived mass limits range from 3.4 to 
8.3 Mj up ; in this interval the more stringent astrome- 
try even indicates the absence of objects from 1.1 to 
0.22 Mj up , i.e. below Saturn mass for periods > 370 d. 

4. Hence no massive planets > 3.7 Mj up exist in orbits 
with periods of 0.75 - 600 d, i.e. at 0.008 - 0.69 AU. 

5. At periods > 600-3000 d (separations > 0.69-1 AU) 
RV derived mass limits range from 8.3 to 22 Mj up . 

6. At the level of our measurement precision the RV data 
are not notably affected by stellar activity. 
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